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ABSTRACT 


The established concept of temperature determination from dual spectral band 
radiometric measurement is analyzed for the case that the measurement device is sensi- 
tive to total photon emittance rather than radiant emittance. A temperature dependent 
ratio of measurements is calculated for two distinct spectral bands using a black body 
distribution. It is shown that, if the bands are nor self contained, then there exists a 
One-to-one correspondence between temperature and ratio. 

A prototype algorithm is proposed and tested which demonstrates the feasibility of 
extracting both the temperature and range by using three distinct spectral band meas- 
urements. The model assumes nonhomogeneous, wavelength-dependent atmospheric 
attenuation. The target and sensor systems each have an arbitrary location in a vertical 
plane relative to the earth. In the computer simulations, LOWTRAN data was used. 
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I. INTRODUCTION 


A. REVIEW OF PRIOR WORK 

In the photon detection process electromagnetic radiation causes three kinds of 
mechanism such as photon effects, thermal effects, and wave interaction effects 
[Ref. 1]. The focus of the remaining discussion will be confined to the first two effects 
which cover the majority of applications involving the detection of thermal radiation. 
One common photon effect is photoconductivity. Photoconductivity is observed in 
semiconductors and characterized by photon energy and energy gap. One common de- 
vice which uses a thermal effect is the bolometer. Bolometers can be made of any ma- 
terial which has a temperature-dependent resistance. This effect is characterized by the 
change of resistance from the heating effect of incident radiation. Therefore, Plank’s 
equation for radiant emittance can only be applied directly to thermal effects. Although 
some devices have a detection mechanism based on the thermal heating effects due to 
radiation, this is not true of semiconductor-based detectors. Semiconductors devices are 
characterized by an energy gap. Only photons with energies greater than this gap can 
produce electron hole pairs which contribute to the detected signal. The excess photon 
energy is dissipated as heat [Ref. 2] and therefore does not contribute to the electrical 
signal generated. Specifically the integrated spectral emittance [Ref. 3], 1.e., photon flux, 
is @ more appropriate signal indicator than the radiant emittance for measurements 
based on semiconductor devices. 

R.B. Johnson and E.E. Branstetter presented a numerical method for the integration 
of Plank’s equation [Ref. 4]. Johnson derived an approximated equation for effective 
radiant sterance [Ref. 5]. Normalization of spectral bandwidth, for convenience in 


analysis, was provided by F.E. Nicodemus [Ref. 6]. Based on Plank’s equation for ra- 


diant emittance, M.]H. Horman presented the technique for temperature determination 
using the ratio of output from two bands [{Ref. 7]. On the other hand, analysis of the 
dual spectral band temperature measurement for devices based on the photon effects has 
not been developed. 

One standard inethod for the determination of the gray-body emissivity and the 
temperature of an object is based on two separate narrowband radiometric power 
measurements. In combination with Plank’s radiation equation for a gray-body, two 
narrowband measurements are sufficient to produce a unique value for both temperature 
and emissivity. For practical reasons, such as improvement of signal-to-noise ratio, 
most instrumentation take measurements over bands which are too wide to be consid- 
ered narrowband. In this more general situation the question of uniqueness of temper- 
ature requires more careful consideration. This issue has been recently analyzed for the 
distribution based on a Plank’s spectral radiant emittance [Ref. 8]. The authors devel- 
oped an analysis based on the evaluation of the ratio (#) of two measurements of the 
integrated power distribution over separate bands. Three cases, nonoverlapping bands, 
partially overlapping bands and totally overlapping bands, were considered. The analyt- 
ical investigation demonstrated that the last case does not always generate a one-to-one 
correspondence between the ratio @ and the temperature. 

For the passive ranging technique, P.J. Ovrevo and R.C. Wood presented the 
method using absorption property of the atmosphere [Ref. 9], and J.R. Jenness, Jr., and 
F.J. Shimukonis proposed the principle of passive ranging using the ratio of the signals 
in narrow spectral bands received [Ref. 10]. 

For optical properties of the atmosphere, R.A. McClatchey and others described 
effects of index of refraction in the atmosphere [Ref. 11: pp. 41]. It was found that 
transmission factors are affected by no more than | - 2 % if refraction is included in the 


calculation. 


B. OVERVIEW OF THE THESIS 

Chapter I is the introduction and covers both a review of relevant literature and a 
thesis overview. In the interest of keeping this thesis fairly self-contained, Chapter I] 
provides a brief technical description of essential background preliminaries. In Chapter 
III, effects due to atmospheric attenuation and detector responsivity are ignored. Using 
dual spectral measurements, it 1s shown that it is possible to uniquely determine a tar- 
get’s temperature if the detector bands are not 100 % overlapping. Specifically, the ratio 
of the dual spectral band measurements 1s shown to exhibit a mathematically predictable 
one-to-one correspondence with the temperature, as long as the bands are not totally 
contained within the other. 

In Chapter IV, a numerical program for calculating a spectral band measurement 
without the assumption of constant atmospheric attenuation and detector responsivity 
is discussed and tested In the process of this development, a novel definition for a 
“responsivity, applicable to devices sensitive to photon emittance, has been introduced. 
If the effects of atmospheric attenuation on the dual spectral band measurement are not 
equivalent, then the ratio scheme for temperature determination, discussed 1n Chapter 
III, is not applicable. It is precisely under these conditions that by using three distinct 
bands both the temperature and the range can be determined. 

In Chapter V, a prototype algorithm is discussed, which demonstrates the feasibility 
of obtaining both the temperature and the range of a target using three distinct spectral 
band measurements. Complications introduced by not assuming constant altitude of the 
target-to-sensor trajectory are addressed. LOWTRAN data was also used here to char- 
acterize atmospheric attenuation. Conclusions for the thesis are presented in Chapter 
VI. In order to not obscure the main points of the thesis, the extensive mathematical 
derivations have been relegated to appendices. Both Appendix A and Appendix B fit this 


description. Appendix C, Appendix D and Appendix E have been created for those 


readers who are interested in an algorithmic description of programs used in the com- 


puter simulation. 


ll. BACKGROUND 


In this chapter, an alternate approach based on photon emittance rather than radi- 
ant emittance is applied. This kind of treatment is appropriate for photovaltaic or 
photoconductive semiconductor devices but not radiometers. First, a few basic as- 
sumptions need to be presented. The discussion starts from the ideal case with the fol- 


lowing preliminary assumptions [Ref. 8: pp. 1256]: 
® opaque object 
e Lambertian surface 
® object in thermal equilibrium 
® spectrally constant emissivity (gray body) 
¢ photon emittance derived from Plank’s radiation law 


¢ unity spectral transmittance of intervening media between the object and measuring 
instrument. 


Plank’s radiation law gives the spectral radiant emittance, which is given by 


fever. 3: pp. 35]: 


2mne(A)he* 


ee. ol scm” spe") 2.1 
[exp(he/AkT) — 1] SCTE (2.1) 


WT) 
where e(/) is the spectral emissivity, h is Plank’s constant, J is the wavelength, c is the 
speed of light, k is Boltzmann’s constant, and T is the object’s temperature. In this sec- 
tion, since the emissivity is assumed to be spectrally constant, e(/) is denoted as «. 

The photon emittance is derived from (2.1) dividing by hAc/A, which is the energy 
associated with one photon [Ref. 3: pp. 38]. Therefore photon emittance is the flux re- 
lated to the number of photons. The resultant photon emittance equation is obtained 


as 


LIEC 


= hotonses-' ecm’ um Lie 
@; A lexp(he[AkT) — 1] P — 


The integral of (2.2) over some band limits gives the photon emittance in the bands. 


The photon emittance within the spectral band limits [4,, 4,] 1s given by 


Ay 
Ody, A, TD) = I O(a, Ta (2.3) 


From the assumptions, emissivity is assumed to be spectrally constant. Therefore 
photon emittance for a gray body is simply expressed as photon emittance for a black 


body multiplied by a constant emissivity and given by 
Q, = €Qpp. (2.4) 


From (2.3), it follows from (2.4) that 


Ay 
Ola dv =e} oad, Da 


= €Opp(42, 41,1), 


(2.5) 


where A, < A, and Q,, (4, T) is the spectral photon emittance for a blackbody [e(4) = 1]. 
It has been shown that the temperature of an object can be determined by taking a 
ratio of measured radiant emittance of an object in each of two spectral bands. 
[Ref. 8: pp. 1256] In Chapter I] and Chapter III, this concept is applied in a simular 
manner except the assumption is made that a ratio is formed from the measured photon 
emittance. Let the two spectral bands be denoted as [/,, 4,] and [/,, 4,] . Let the ratio 


for photon emittance be defined as 


O(A,, 4, T) 


GU dein eo) 


BA, Aas A3, Aa, T) = 
Since from the assumption, emissivity is uniform over the spectral band, emissivity 1s 
cancelled out, and the ratio becomes the ratio of the photon emittances associated with 


a black body, as shown by 


OnplAd, Ay, T) 


BA, 49,43, 44, =o. 
Ay Zoey oon T) Ope yew ees @) 


(2.7) 


This implies that @, is a function of the object’s temperature T only and the ratio can 
indicate the object temperature. When the dual spectral-band photon emittance method 
is used, it 1s critical to know that the ratio @, and the temperature have a one-to-one 
correspondence. For this purposes, three general cases of spectral bands which cover 


all possibilities are considered. The possible cases are [Ref 8: pp. 1256] 
1. non-overlapping bands (Figure 1a), 
2. partially overlapping bands and common band limits (Figure 1b,c), and 


3. totally contained band (Figure 1d). 


These three cases are presented graphically in Figure 1. In the following two chapters 
the validity of temperature determination based on dual spectral photon emittance for 


the ideal case is discussed. 


a. nonoverlapping case 


1 2 3 4 
b. partially overlapping case 


1 3 2 4, 


Cc. one common band limit case 


1,3 2 4 
d. totally contained case 


Figure 1. graphical representation of possible band combination: three cases are 
nonoverlapping, partially overlapping, and totally contained. 


Il, DUAL SPECTRAL BAND MEASUREMENTS FOR TEMPERATURE 
DETERMINATION 


In this section the ratio of photon emittance defined in (2.6) is considered for the 
various cases discussed at the end of Chapter II. In order to validate the dual spectral 
bands method, it is required to show one-to-one correspondence between the ratio and 
the temperature of the object. For this purpose, the derivative of the ratio of band lim- 
ited photon fluxes, @, in (2.7), with respect to temperature is performed. If the ratio of 
photon emittance @, is proved to be either strictly increasing or strictly decreasing 
function with respect to temperature, 1.e., either 02,/0T > 0 or 0&,/0T <0 for all T, then 
the ratio of photon emittances and the temperature exhibit a one-to-one correspond- 
ence. It follows that this approach is potentially applicable to the temperature determi- 
nation of grey bodies. 

To simplify the integral of (2.2) , the variables of integration in photon emittance 


QO(4,, 4,,1) are replaced by the dimensionless parameters [Ref. 8: pp. 1257] 





Pe ue 
and 
= Li 32 
4 Fae T 9 ( v a) 


where TJ, is an arbitrary reference temperature. By prior convention [Ref. 8: pp. 1257], 


it is chosen to be 


T, = (3.2b) 





where /,, as introduced in Chapter I], is the highest wavelength in spectral band 1. Using 


(3.2b) as a reference temperature in (3.1), the spectral band limits of w are given by: 


he 
AKT, 





x; = (sea) 
where /; is one of the spectral band limits [Ref. 8: pp. 1257]. After substitutions into 


(2.5), the new representation for the photon emittance is given by 


Q(2), 4), T) = 2nek? eh | (Wy, (3.4) 


aX; 


where [Ref. 8: pp. 1257] 


y? 
g(y) = i) a (3.5) 
The ratio of photon emittance is obtained as 
[sway 
Z,= = ea (3.6) 
| g(y)ay 


and the derivative of the ratio of photon emittance signals (3.6) with respect to temper- 


ature comes out to be 


[rae [ewnav - | "sonay [nna 
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(3.7) 
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where 


3 
fl) = exp(y)¥ 


_ ; 3.8 
[exp(w) — 1} - 


(The derivation of (3.7) is in Appendix A.) To simplify the expression of (3.7), the fol- 
lowing operator is defined [Ref. 8: pp. 1257]. 


ax, aX, ax, aX, 
[a.bI, x», 22,260) =] away] “oway -[owrav | awry (3.9) 


Then the ratio becomes 





(Agid, Xx; X3; X4; a) 


| “ue 


aX3 


l 
aT 1 (3.10) 


Since the denominator of (3.10) is always positive, the sign of 0#,/0T is determined by 
the sign of the operator [f,g]. It is necessary to show that the sign of [f,g], which is the 
sign of (3.10), does not change with a, or equivalently with temperature. The temper- 
ature dependence of [f,g] appears only in a, which can be seen in the limits of integration. 


The integrands do not have an explicit dependence on temperature. 


A. NON-OVERLAPPING CASE 
In order to prove that the sign of [f,g] is fixed over the domain of a, several proper- 


ties of ftw) and g(W) are required. They are [Ref. 8: pp. 1257] 
1. /(0) = g(0) =0 
2. 0<g(W) </ly) for y >0 
3. lim,_..W) = lim,_...g(W) = 0 
4. For B>1,f(W) — Bge(W) has exactly one root on the interval w é (0, oo). 


The discussion to follow will clarify the significance of these requirements. Graphical 
representation is useful to clarify these properties and it is shown in Figure 2. This figure 
is especially useful for visualization of the fourth property. For values of # > 1, it is ob- 
vious that f(W) crosses Bg(w) exactly once. A more formal discussion of the four prop- 
erties is provided in the next paragraph. 

The first property is proven by using l’Hopital’s rule. In order to prove the second 


property, define 


PCW) = ae (3.11a) 


which by direct substitution from (3.5) and (3.8) is given by 


P(p) = (3.118) 


a 
= eye 


It then can be shown that P(w)>1 for ¥>0. This follows from the derivative of P(w), 


which 1s: 


dP(y) __1—exp(— W)— vexp(-W) 


Jah? 
ay (1 = exp(— v))" = 
After expanding (3.12a), it is clear that : 
GPW) (3.12) 


dy 1—exp(—¥) (1 —exp(—)/ 


which follows from the observation that the second term in the expansion is smaller in 
magnitude. Also, by direct application of I’'Hopital’s rule on (3.11b), lim,_,-P(W) = 1. 
These two results imply that P(W) > 1. Thus /(W) > g(W). Consequently, after noting that 


both /(W), as given in (3.8), and g(w), as given in (3.5), are strictly positive, property two 


follows. The third property is simply proven by limit evaluation. For the fourth prop- 


erty, it is necessary to evaluate f(w) — Bg(W) given by 


___exp(yye v 


which can be expressed as 


y w exp(y) 
Sy) — Ba) > Seen a p} (3.130) 


By noting (3.11b), it follows that 
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f(b) — Bg(v) = 7 LP) — BI. (3.13c) 


W 
exp(y) — 


Since lim,_,-Q(w)=1 and B is constrained to be greater than one, the quantity 


Sl) — Be() satisfies, after reference to (3.11a), the inequality 
Sv) — Bgl) <0 (small y) (3.14a) 


for sufficiently small yw. Also since lim,_.,.P() = oo , the inequality has to reverse for 


large enough w. In other words 
Sv) — Bg(v) >9 (large ¥). (3.145) 


Since, from (3.12a), P(W) is a strictly increasing function, it is correct to conclude the 
inequality reversal can only occur once. Property 4 is therefore proven. 
For the non-overlapping case represented on Figure la, there are only two sets of 


relations which can define the spectral bands (3.3): 


Ns nx, (37552) 


and 


x3 < X4 < l< X4- (3.155) 


Relation (3.15a) is considered first. Motivated by the fourth property previously stated, 
let 
Slax) 


Be) eae (3.16) 


The quantities f(w) and B(a)g(w) must cross each other only once at the value of wy de- 
noted as x, This is represented in Figure 3. The geometrical aspects of figures dictates 


the following relations for all values of a: 


[Pavan < | aeewray (3.172) 
and 

aXy aXy 

[“rnrav > [“remanav. (3.178) 


Substitution (3.17a) and (3.17b) into (3.9) shows the operator [f,g] satisfies: 


aX? aXg aX? aXy 
(felll, x23, 20 0) < |” away | 'ewiay — [awa] slode(v)av = 0 


(3.18) 


for all the values of « (3.2a). According to the brief description following (3.10), this 
proves that 0#,/0T <0 for all the values of T and therefore &, is a strictly decreasing 
function of temperature. For the second relation of non-overlapping bands (3.15b), a 
similar analysis would show that 0@,/oT > 0 and @, 1s a strictly increasing function of 


temperature. 


B. PARTIALLY OVERLAPPING AND COMMON BAND LIMIT CASE 
In this section, the case which is partially overlapping and the common band limit 


case are considered. Two sets of relations are considered in this case and they are 


l<x3<x,< x, orl <x3< x, 5x4 (3.19a) 


x35 1<x4< x, o1rx3< 1 <x sx, (3.195) 


For this discussion, the graphical method previously introduced is useful [Ref. 8: pp. 
1258]. Figure 4 illustrates the first set of the relations (3.19a) for partial overlap. With 


reference to Figure 4, six sets of integrals are defined as follows: 


s= [nwa 4.= [ewe (3.204, 

=| Andy A= | enay (3.2064) 
aXy aXy 

s= [Away =| eWav. (3.20e/) 


After inserting this formalism into (3.9), the bracketed operator [f,g] becomes 
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ell, x2, 3, %q, &) = (S, + S)(A2 + Aj) — (A; + A2)(S2 + S3)- (3.21a) 


After expansion and recombination of (3.21a), it follows 


Ag l(l, x2, 43, %4, a) = [A3(S, + S2) — S3(4, + A2)] + (S142 - A,S). (3.215) 


This can be retranslated using (3.20) into two terms which follow 


[fgl(1, x2, x3, 4, %) = [ig l(1, x2, x2, %4, a) + (Ael(l, X35 X3, Xp, O). (3.21c) 


Case (3.19a) will now be applied on expression (3.2Ic). If 1<x,<x,<.x,, then both 
terms in (3.2Ic) are of the form | < x, = x; <x, which satisfies the nonoverlap condition 
(3.15a). Since, as indicated in the previous section, the particular nonoverlap condition 
(3.15a) generates negative values for [f,g], the expression (3.21c) must be negative for the 
stated inequality. 

For common band limit case, which is 1 = x, < x, < x, or 1 <x,;< x, =x,, one of the 
terms 1n (3.21c) becomes zero. However, the other term remains negative, and therefore 
[f,g] remains negative. This discussion implies that the ratio @, 1s a strictly decreasing 
function of temperature. Similarly, for the second relation of this case (3.19b), the ratio 


#, iS a Strictly increasing function of temperature. 


C. TOTALLY CONTAINED SPECTRAL-BAND CASE 

This section discusses case 3, where one interval is entirely contained within the 
other. By the use of asymptotic approximations [Ref. 8: pp. 1258], 1t 1s possible to show 
that the ratio 2, may have the same value for two or more temperatures. Therefore, the 
measurement to determine the temperature is not unique. Two examples, which are the 
case where « is large and the case where a is small, are discussed. 

As previously discussed, the sign of 0@,/0T is determined by the sign of [f,g]. There 


are two sets of relations which are 1 <x,<x,< x, and x,;< 1<x,<.x,. The first set of 


relations is arbitrarily chosen for discussion since the approach is similar for the other 
relation. For large or small value of a asymptotic methods are useful to determine the 
sign of [f,g]. (The derivation using asymptotic approximation is in Appendix B.) For the 


large value of a (i.e., small temperature ), [f,g] is approximated as 
fg]xoxs (1 — x3) expL — a(x, + 1) J. (3.22) 


Since x, is greater than one, [f,g] becomes negative. For small values of  (..e., large 


temperature ), [f,g] 1s approximated as: 


[hele a0 x3)(x} — 1) — (x9 -1)(g — 34) J. (3.23) 


Given x, and x,, it is always possible to find a sufficiently large x, such that (3.23) 1s 
positive. Combining the two results for small a and large a, the sign of 06@,/0T is not 
fixed over @ ({i.e., temperature ). This means that the ratio 2, does not have a one-to-one 
correspondence with temperature. Therefore, for the case when one spectral band 1s 
completely contained within the other, the temperature determination method based on 


the dual spectral band measurement of photon emittance is not, in general, valid. 


D. RATIO CALCULATION FOR PHOTON EMITTANCE 

In order to generate the output ratio calculation, the numerical method for inte- 
gration known as Simpson’s rule [Ref. 12: pp. 95] was used. The flow chart for this 
calculation is shown in Appendix C. The examples were demonstrated using the condi- 
tions in [Ref. 8: pp. 1258-1259]. One case has a band combination of [2.5 um, 5.5 pn] 
and [3.9 um, 4.1 um] shown in Figure 5. The second case has a band combination of 
[10.0 wm, 14.4 pm] and [11.8 um, 12.2 ym] shown in Figure 6. In these figures, the nor- 


malized ratio denoted as @, was used to enhance the results. Normalization is done by 
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2, = 3.24 
i Dor = heel ( ) 


where Z,,,, 1s the maximum value of &, and @,,,, is the minimum value of &, Both cases 
are slightly different from the results in[Ref. 8: pp 1259], but the general features are the 
same. One ratio value may correspond to two or more temperature values as predicted 
in Section C. These figures indicate that the minima of the ratio occur in different po- 


sitions depending on the bands selected. 


E. DISCUSSION OF COMPUTER RESULTS 

Computer simulations were performed for the three cases discussed in this chapter. 
The first example is the nonoverlapping case and the band pair is chosen to be [3.0 um, 
4.0 pm] and [4.0 um, 5.0 pm]. Figure 7 shows this example and it exhibits the decreasing 
feature as expected in Section A. The second example is the partially overlapping case 
and the band pair is chosen to be [3.0 um, 4.5 wm] and [3.5 um, 5.0 um]. Figure 8 shows 
this case and the results agree with the discussion in Section B. The third example is the 
case of totally contained band within the other and the band pair is [3.0 um, 5.0 um] and 
[3.5 wm, 4.5 um]. Figure 9 represents this case and the results show that one-to-one 
correspondence does not exist in this case. It can be seen from Figure 9 that a ratio value 


of 2.0 corresponds to both 460 degrees Kelvin and 1400 degrees. 





Figure 2. Plot of f(y) and Be(y) for B= 1,2,...,5 





Figure 3. | Nonoverlapping case.: The constant f(a) is chosen such that 


B(a)g(ax,) = flax). 
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Figure 4. Partially overlapping case or common band limit case: Each region 1s 
defined in the text. 
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Figure 5. Example of totally contained case (normalized ratio): Bands are [2.5 
um, 5.5 pm] and [3.9 um, 4.1 pm]. 
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Figure 6. Example of totally contained case (normalized ratio): Bands are [10 
um, 14.4 wm] and [11.8 wm, 12.2 um. 
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Figure 7. Example of the non-overlapping case: Bands are [4.0 um, 5.0 um] and 
[3.0 um, 4.0 um]. 
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Figure 8. Example of the partially overlapping case: Bands are [3.5 ym, 5.0 um] 
and [3.0 um, 4.5 pm]. 
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Figure 9. Example of the totally contained case: Bands are [3.0 um, 5.0 um] and 
[3.5 um, 4.5 um]. 
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IV. A GENERAL MODEL FOR TEMPERATURE DETERMINATION VIA 
DUAL SPECTRAL BAND MEASUREMENTS 


In Chapter II, ideal photon emittance temperature determination was discussed. 
Since this discussion was dealing with the ratio without considering measurement of the 
photon emittance as a system, it was not always practical. Therefore, some other factors 
should be taken into account to realize this method. They are emissivity of the object, 
atmospheric transmittance, detector responsivity, and so on. In the previous discussion, 
those parameters were assumed to be spectrally constant or unity. Strictly speaking, 
they depend on the wavelength. Also after considering these factors as part of the sys- 
tem, the resultant ratio becomes the ratio of the output voltages from the detectors of 
two bands instead of the ratio of photon emittance. A brief discussion of these param- 


eters follow. 


A. EMISSIVITY 

The emissivity dictates how much thermal energy is radiated from an object. Its 
value lies between zero and unity depending upon the material characteristics of thermal 
source. It is denoted as e(2). Although emissivity for some materials has temperature 
dependency, it is assumed to be independent of temperature to make discussion simpli- 
fied. Materials are divided into three types based on the value of emissivity. They are 


[Ref. 3: pp. 40] 
1. A black body, e(A) = I. 
2. A gray body, e(/A) = constant < 1. 


3. A selective absorber, e(/) varies with wavelength. 


The spectrally constant emissivity cases were covered in the previous Chapters II and 


II]. Therefore, wavelength-dependent emissivity is the primary focus of discussion in this 
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chapter. In addition, since emissivity depends on the types of materials, it is not possible 
to get reliable analytical expressions for all selective absorbers. Usually emissivity values 


are obtained by measurement. 


B. ATMOSPHERIC TRANSMITTANCE 

The atmosphere is known to behave as a wavelength-selective absorber. This imn- 
plies that the transmission of signals through the atmosphere becomes wavelength de- 
pendent. Specifically, the atmosphere affects the transmission of infrared radiation. The 
main compositions of the atmosphere are nitrogen, carbon dioxide, oxygen, and water 
vapor. Among these, molecules which have nonsymmetric molecular structure are pri- 
mary responsible for absorption. Thus carbon dioxide and water vapor are the dominant 
sources for absorption. In addition some other small constituents of the atmosphere, 
which have non-symmetric structure, also contribute absorption. 

Atmospheric scattering is the other important phenomenon which affects atmo- 
spheric transmission. The strongest atmospheric scattering occurs in the case when the 
size of aerosols and the wavelength of the signal are almost the same. In this particular 
case the scattering coefficient depends highly on a wavelength of the signal. The Mie 
theory describes this phenomenon. 

The atmospheric transmittance is determined by atmospheric absorption coefficient 


and atmospheric scattering coefficient. It is given by 


t(A, x) = exp[— (a(A) + y(4))z] (4.1) 


where a(/) is atmospheric absorption coefficient, y(A) 1s atmospheric scattering coeffi- 
cient, and x is the distance between the object and the detector. An example 
transmittance data by the NPS LOWTRAN6 program [Ref. 13] at 20 meters in height, 
0 degree elevation for 1 km and 5 km distances under 1962 U.S. standard atmospheric 


conditions are shown in Figure 10. The importance here is that the transmittance de- 
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pends not only the wavelength but also distance and atmospheric conditions. Moreover, 
atmospheric conditions vary with height and thus the transmission depends on the path 
between the object and the detector. Unless the path and atmospheric conditions are 
Known, it 1s impossible to apply the ideal ratio calculation. Therefore, the transmission 
path and atmospheric conditions have to be provided in order to obtain the atmospheric 
transmittance. A recent report [Ref. I1: pp. 41] on transmittance calculation for the at- 


mosphere has concluded that refraction effects can essentially be ignored. 


C. DETECTOR RESPONSIVITY 
The detector responsivity 1s generally obtained by dividing the measured output 


signal voltage by the optical power falling on the detector. And it is defined as 


Ve 


TA: 





(4.2) 


where R is the responsivity [V/W], V,,, 1s the RMS output voltage, H is the RMS 
irradiance [JV/cr?] at the detector location, and A, 1s the detector area [cm]. However, 
this quantity is not the appropriate parameter for the development of the photon 
emittance case. A new proposed responsivity R’ is defined by modifying (4.2) in order 
to adjust it for the photon emittance calculation. It is given by dividing the measured 


output voltage signal by the number of photons into the detector. The proposed R’ is 


defined as 


R' — Vas (4 3) 
OAg 


where R’ is the photon responsivity [V+ photons-'+ s], and Q is the photon emittance 
[photonses'+cm]. It follows from definitions (4.2) and (4.3) that the photon 


responsivity R’ is simply obtained by multiplying responsivity R by he/A, 
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Ri = RE, (4.4) 
as shown above. For photon detectors sensitive to the number of photons, R’ should 
be fairly constant out to the cut off wavelength. 

Unfortunately, detector responsivity (4.2) is not usually provided among detector 


specifications, but the specific detectivity D* is given instead. The D* is the normalized 


quantity of detectivity defined by [Ref. 3: pp. 270]: 
D* = DJ AGB (4.5a) 


where B is the electrical system bandwidth in heltz, A, is the area of the detector, and 
D is the detector detectivity. This quantity is often expressed in terms of the reciprocal 
of the noise equivalent power (NEP). Therefore, the D* is also given by [Ref. 3: pp. 
270] 


JAB 
x 
De (4.56) 





The noise equivalent power is the radiant flux required to obtain an output signal equal 


to the noise of the detector. It is given by 





(4.6) 





where V, is RMS signal voltage and Vy is RMS noise voltage. From (4.6), the specific 


detectivity becomes 


Dt = (4.7) 
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The photon responsivity is obtained by combining (4.3), (4.4), and (4.7), and it is 


D*V 
Re le (4.8) 


JAgB 4 





This equation shows the specific detectivity is closely related to the photon responsivity. 
From the expected wavelength dependencies in (4.8), it follows that the quantity D* 
should increase linearly with wavelength out to the cutoff wavelength [Ref. 3: pp. 294]. 
Since photon detectors respond to the number of photons absorbed in the detector up 
to cutoff wavelength (usually dictated by the band-gap), the modified responsivity of 


photon detectors is ideally independent of wavelength. 


D. COMPUTER SIMULATION OF TEMPERATURE DETERMINATION 
Based on above discussion, the ratio of photon emittance for black body is replaced 


by the ratio of the output voltages at the detectors, and it is given by 


[Peay oR awa 
ay =~ ——_________ (4.9) 
[Mea oR eWay 


a 
where é(/) is the emissivity, t(/, x) is atmospheric transmittance, and R’‘(/) is the photon 
detector responsivity. The subscript on the photon detector responsivity denotes the 
spectral band used in the measurement. In this equation, assuming that the wavelength 
parameters are spectrally constant and two identical detectors are used, the ratio be- 
comes identical to the ratio of photon emittance for black body (3.6). In order to gen- 
erate the output ratio calculation under more general conditions, the numerical method 
for integration known as Simpson’s rule [Ref. 12: pp. 95] is used. The algorithm for this 


calculation is shown in Appendix D. 


ot 


1. Discussion of the computer simulations 

As described in this chapter, the actual ratio for the photon detector is the ratio 
of the output voltages. In the rest of this section, the simulation for the ratio calculation 
is discussed. Although the photon detector responsivity has wavelength dependence, as 
discussed in Section C, it is not a very sensitive factor up to the detector cut-off. 
Therefore, this is ignored in this discussion, and emissivity and atmospheric 
transmittance are considered. Since emissivity depends on the material of the thermal 
source, many kinds of emissivity values are possible. For this simulation, an arbitrary 
created emissivity curve shown in Figure 11 was used. For atmospheric transmittance, 
the atmospheric conditions of 1962 U.S. standard atmosphere were selected. The data, 
shown in Figure 10, are obtained from NPS LOWTRANG6 program. For the purpose of 
comparison with the results for ideal case, the same band pairs are selected as ones in 
Chapter III. The result for overlapping case 1s shown in Figure 12, for partially over- 
lapping case in Figure 13, and totally contained case in Figure 14. The curves generated 
in these figures correspond to | km and 5 km. With the exception of Figure 14, there 
is a One-to-one correspondence between measured ratio and temperature for each range. 
In Figure 14, there is a slight positive shape after T=1300 K which violates the 
uniqueness condition. Although the values are shifted slightly, the shapes of the curves 
are essentially the same as ones for ideal ratio case shown in Chapter III. The general 
features of each case are the same. However, the results demonstrate that the curves 
depend on distance. This is quite inconvenient. The application of two-band temperature 
determination is then impossible without range information. Also, to lesser extent, the 
application is limited since the emissivity is material dependent. If the source for thermal 
radiation is not a gray body, then the emissivity of the source needs to be known 1n ad- 


vance. 
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Figure 10. Atmospheric transmittance for 1962 U.S. standard atmosphere: The 
path lengths are taken to be 1 km and 5 km in the horizontal direction. 
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Figure 11. Example of the emissivity (arbitrary chosen): This data does not cor- 
respond to the specific material. 
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Figure 12. The ratio of output voltages for nonoverlapping case: Bands are (5.0 
um, 4.0 um) and (4.0 um, 3.0 um). The distances are 1 km and 5 km. 
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Figure 13. The ratio of output voltages for partially overlapping case: Bands are 
(5.0 um, 3.5 wm) and (4.5 um, 3.0 um). The distances are | km and 5 
km. 
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Figure 14. The ratio of output voltages for totally contained band case: Bands are 
(5.0 um, 3.0 wm) and (4.5 um, 3.5 wm). The distances are | km and 5 
km. 
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V. THERMAL RANGING TECHNIQUE 


From the previous chapter, the ratio of output voltages from dual spectral bands is 
not the same as the ratio of black body photon emittances under realistic and nonideal 
conditions. To solve and apply the method under these conditions can be quite difficult. 
On the other hand, one significant idea is hidden in this problem. 

The concept presented in this chapter is created from the difficulty, in the temper- 
ature determination method, produced by nonconstant atmospheric transmittance. 
Among assumptions discussed in Chapter IV, atmospheric transmittance breaks the va- 
lidity of temperature determination procedure because it depends on wavelength, the 
distance between an object and the detector, and atmospheric conditions. Although 
emissisivity and detector responsivity have to be taken into account, under some rea- 
sonable assumptions they can be neglected. It will be shown that the dependency of at- 
mospheric attenuation on the distance makes the concept of thermal ranging feasible. 

The use of infrared thermal signals has the advantage that active interrogation 1s not 
required since the object radiates energy or photons. This passive property of infrared 
detection can be significantly important in military situations. However, the passive de- 
tection of the target does not automatically provide range data and this is a disadvantage 
of thermal infrared detection systems. One interactive (nonautomatic) passive system 
involves the use of a trained human observer. The training allows the observer to rec- 
ognize a target and by noting the size of the target on the screen make an experienced 
estimate of the range. Currently ranging is done by active sensor systems such as radar 
or laser range finder. These active systems are capable of being sensed by the target. 


Although active ranging systems tend to provides precise distance information, they ex- 
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pose the instrument generating the signal. Therefore, the development of a passive 
ranging system would provide significant military advantage. 

The Electronic Support Measure (ESM) equipment in electronic warfare systems is 
the passive sensor in the microwave region. This gives target identification, direction of 
the target and also a very rough range information. For this method, the power trans- 
mitted by a radar of a target, the antenna gain, and attenuation factor must be obtained 
in advance. It seems that the same concept can be applied to infrared systems only if 
the temperature of a target, surface area, emissivity, and atmospheric transmittance are 
known. Therefore, this method 1s not practical. 

The following discussion proposes a new concept of passive ranging by infrared 
systems. It will be shown that this concept 1s a feasible solution to the problem of target 
ranging using only automatic passive methods. The following discussion consists of 


maree Steps: 


1. Step 1: Atmospheric attenuation coefficient is assumed to be constant. Temper- 
ature of the target is known. 


2. Step 2: Atmospheric attenuation coefficient is assumed to be constant. Temper- 
ature of the target is unknown. 


3. Step 3: Atmospheric attenuation coefficient is not constant. Temperature of the 
target 1s unknown. 


The concept of ranging shown here is based on the distance dependency of atmospheric 
transmittance. Step | and step 2 are intermediate procedures and they do not allow 
changing atmospheric conditions due to, for example, an elevated transmission path. 
In step 3, these prohibitions are rescinded and the ranging model has greater applica- 
bility. 
A. RANGING FOR KNOWN TEMPERATURE CONDITION WITH 
HOMOGENEOUS ATMOSPHERIC CONDITIONS 

In the analysis of this section the emissivity of the target and photon detector 


responsivity are assumed to be spectrally constant. In addition, it is assumed that at- 


yy 


mospheric conditions are not changing along a line-of-sight trajectory between target 
and passive sensor system. Therefore, atmospheric attenuation coefficients will not de- 
pend on range but still exhibit wavelength dependence. By appropriate band selection, 
i.e., from a flat portion of the transmission curve [LOWTRAN], the wavelength de- 
pendence of the transmission can be approximately neglected with in each of the selected 


bands. Therefore, atmospheric transmittance for the i-th band is given by 


(x) = exp( — 072) (5.1) 


where oa, is the atmospheric attenuation coefficient within the i-th band, and ~x is the 


transmission path length. Applying the assumptions stated to the ratio (4.9) leads to 





gM = R! (5.2a) 


After substitution of (5.1), the ratio can be written as 
AM = exp[— (0; — 0,)x]R’ (5.2b) 


where &” is the ratio of the measured output voltage at the detector, #’ is the ideal ratio 
of photon emittance, and o, are the atmospheric attenuation coefficient for the measured 
bands. In (5.2b), the ratio of the output voltages is obtained by direct measurement. The 
atmospheric attenuation coefficient is calculated from the atmospheric conditions, and 
the ratio of photon emittance is determined by the temperature given. Therefore, the 
transmission path length is the only unknown parameter. From (5.2b), the transmission 


path distance 1s obtained as 


al ee (35:3) 
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The use of (5.3) makes passive thermal! ranging feasible. However, it is applicable to very 
limited conditions. In the battle field, the targets are usually unknown, and consequently 
the temperature of the target is hard to obtain. This means that the system should be 
capable of determining both range and temperature to be practical. The new proposed 


concept is described in the following sections. 


B. RANGING FOR UNKNOWN TEMPERATURE WITH HOMOGENEOUS 
ATMOSPHERIC CONDITIONS 

In this section a concept to determine the range and the temperature of the target 
is discussed. As in Section A, the atmospheric conditions are assumed to be homogene- 
ous in space. Mathematically, two nondegenerate equations can be used to solve for two 
unknowns. In this discussion, the two unknowns are range and temperature of a target. 
To get two equations, two distinct sets of ratio calculations are required. For that pur- 
pose, no less than three sets of band measurements ( 1.e., band 1, band 2, and band 3 ) 
are required to produce two distinct equations such as (5.3). All the assumptions de- 
scribed in the previous section are also applied in this section except that the temper- 
ature of a target is unknown. The mathematical derivation to get the range equations 
is exactly the same as in the previous section. The two range equations for three sets of 


band measurements are 





M 
] Zaye. 

x= =e In (5.4) 
GM 


where the subscripts serve as an index to identify spectral bands. Combining (5.4) and 


(5.5), one equation to solve the temperature 1s obtained as 
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(5.6) 
However, this equation is impossible to solve analytically. A numerical method can be 
used to solve this equation. After getting the temperature, the range is easily obtained 
by using (5.4) or (5.5). Since the assumption that the atmospheric conditions are ho- 
mogeneous in space is only approximately valid for very short ranges, this method is not 
practical. The following section solves the difficulty of the determination of the range 


without assuming homogeneous atmospheric conditions. 


C. RANGING FOR UNKNOWN TEMPERATURE WITH INHOMOGENEOUS 
ATMOSPHERIC CONDITIONS 
For the following discussion, several equations are recalled for convenience. For 


four bands, the ideal ratios calculated from photon emittance are defined as: 


O(4,, 49, T) 
gi =, 5.7a 
OC —_ 
and 
Os, 46; T) 
R=, 5.7b 
= OC, dg D = 
where the total photon emittance is given by 
A, 
OA; Ai, T) =| OA, T)d2, (5.8) 
A; 


and spectral photon emittance Q, can be expressed in terms of the gray body distribution 
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2rEC 


eee ee 5. 
AI exp(he/kT) — 1] ae 


0, 
With two band measurements, the temperature can be determined. With three band 
measurements, it should be possible to determine both temperature and range. As pre- 
viously discussed, range information is available as long as the effective atmospheric at- 
tenuations in each band are distinct. The problem posed by obtaining range information 
for a thermal source at the same elevation as the receiver (homogeneous atmospheric 
conditions) is not as complex as that associated with arbitrary elevation. Elevation 
changes in the trajectory of the radiation usually require consideration due to the cur- 
vature of the earth. In the general problem to be solved, both the receiver and the ther- 
mal source take arbitrary locations in the plane defined by the vertical line passing 
through the receiver and thermal source. This 1s represented in Figure 15. 
In the algorithm to be discussed, three arbitrary but distinct radiation bands are se- 
lected. Initially disregarding the effects due to the atmospheric attenuation, calculated 


ideal ratios of signals are given by 


Ideal signal Band 1 


Ideal signal Band 2 © 102) 


BT) = 


Ideal signal Band 3 


Ideal signal Band 1 © ee) 


R3(T) = 

Figure 16 shows the temperature dependence of these curves for bands defined by: 
(A,, 45) = (12.0 wm, 10.8 pm) 
(A3, 24) = (10.8 ym, 9.6 pn) 
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(As, Ag) = (9.6 wm, 8.2 pm) 


The atmospheric attenuation in the form of transmittance is provided from the NPS 
LOWTRANG program [Ref. 13] using data corresponding to a vertical path. Atten- 


uation coefficients are derived from 
| 
oj(4)=—- 7 In t,(A) [nepers/m] (5.11) 


where 1= 1,2,3 indicate the wavelength band and / 1s the index of the specific layer. The 
quantity 7, is the length of the trajectory in the j-th atmospheric layer. See Figure 17. 
Each r, has the same length if the line-of-sight angle 1s 90 degree. The attenuation coef- 
ficient o,, is obtained by taking the 90 degree angle. The standard band average for 1, 
is given by 


A 


upper 


«(A)dd 


Aigner 
aS ==: (5.12) 


With reference to Figure 18, the effective approximate transmission factor t,, for the i-th 


band, is given by: 


Kmax 


t, = exp( — eae (5.13) 
k=) 


where &,,,, 1s the number of layers between source and receiver, and m 1s the layer num- 
ber of the receiver. Here, r, is the length of the trajectory in the layer k+m-—1. The 
approximation involved in (5.12) will be discussed at the end of this section. Expression 


(5.13) can be recast into the form: 


t, = exp( — 9,R) (5.14) 


where 


Kmnax 
Fik+m—=1)'k 
0, = =! ? (a)5) 
and when the range R must satisfy: 
Krnax 
R=) 7% (5.16) 
k=] 


With the exception of the first and final laver, the value for each r, 1s obtained from the 


rule 


i—] 
n=R- > (5.17) 
j=1 


as shown in Figure 19. Each R, is calculated from the Cosine Law relation derived from 


the geometry of Figure 20 
(H, + Re)’ = Rp + (hy + Rey — 2Rihg + Re) cos(6 + oy (5.18) 


or equivalently 
(H, + Rey = Rj) + (hy + Re) + 2Rihg + Rz) sin 8, (5.18) 
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where @ is the line of sight angle between receiver and source measured with respect to 
the horizontal. R, can be obtained from (5.18) using the quadratic formula with the 


correctly chosen ’+’ sign. R, is also obtained from (5.18) after taking 
Hy = {int -) + 1A 5.19 
1 = Lin 52) + dr (5.19) 


where the int[_] results in the integer part. This follows from examination of the geom- 


etry in Figure 15. For R,,...,R,, R, values for H, follow 
H, = H,+(k—-1)Ar (5.20) 
which are used in (5.18) to calculate subsequent R,. For clarity, it is noted that: 


R =r, (5.21a) 


R, =n, + 1, ace Aa (5.278) 


until kK =k,,,, as defined by (5.16). Equations (5.17) combined with (5.20) are solved it- 


eratively to generate the values for r, except for r,,. Lastly, r,,., is calculated from: 
Nemax = Vinal = Kx, Siew =I (5.22) 


Once r, are determined for all i up to &,,,, , the effective transmission for each band (5.14) 
can be calculated from (5.15) and (5.16). The measured ratios, i.e., including the effects 


of the atmosphere, are related to ideal ratios (5.7) according to: 
M r fy 
By =Ay2 ey (5.23a) 
and 
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M 1 %3 
R3, = B23) CAL (5.236) 


in agreement with (5.2a). After substitution from (5.14), (5.23) can be written as 
M I 
By7 = Zy2 exp[— (0, — 42)R] (5.24a) 
and 
M I 
R3; = A; exp [— (03 — 4,)RI, (5.245) 


where R is the distance between thermal source and receiver. Figure 21 shows the 


measured ratios calculated from (5.24) for conditions: 
© o,=0.0001494 [nepers/m] 
® o,=0.0001027 [nepers/m] 
® «o,=0.0004622 [nepers/m] 
°° R=100 [ml 


Solving for R from (5.24), it follows that 





g™ 
R'(N = In —“~ , (5.25a) 


and 





RC 1 Ot (5.256) 


where the subscripts are employed to denote spectral bands. In the algorithm, it will be 
assumed that the correct measured ratios, consistent with (5.24), are provided by exper- 
imental testing. Ideal ratio curves can always be calculated according to the methods 


(5.10) previously discussed. Constants o,, o,, 4, can be approximately calculated from the 
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NPS LOWTRAN6 data. By checking for the crossing point between the two curves 
(5.25), both the temperature, T, and the range, R, are uniquely determined as shown in 
Figure 22. 

The algorithm sketched out is not yet complete. The values of ¢,, 0, and oa, are de- 
pendent on the range, R, since the weights for evaluating the coefficients a; (5.15) depend 
on all the values for 7,. In order to calculate the coefficients a, (5.15), an initial guess is 
made for a range candidate, R,,, within limits dictated by the domain of LOWTRAN 
data sampled. This value for R,, is used in the calculation of o,,¢,, and o,. The search 
for a numerical crossing point for curves (5.25) produces a new value for R which can 
be denoted as R,,. The value for Rg, generated by this first iteration is expected to be 
closer to the actual range but with an unacceptable level of error. The process is re- 
peated by iterating with the R,,. The next estimate generated by the algorithm, Ro, 1s 


tested against the previous one, R,,. In general, the procedure is stopped when 


| ae 7 Rem-1 
Gr 





<6. (5.26) 


where 6 is the predetermined convergence parameter and the index m specifies the 
number of iteration cycles. In Appendix E, the algorithm for the process 1s presented. 
For the purpose of checking the algorithm, the range candidate, Rg, is selected as 
shown in Figure 23. Actual ranges are chosen from 10 meters to 700 meters and Rq's 
are 700 meters for up to 350 meters of actual ranges and 10 meters for longer than 350 
meters. This choice appears to correspond to worst-case guessing. The results for error 
calculation are obtained for 6= 0 degrees, 30 degrees, and 45 degrees at a detector 
height of 5 meters. The actual temperature and convergence parameter 6 are chosen to 
be 353.54 K and 0.001, respectively. The error of range determination is shown in Figure 


24, and the error of temperature determination is shown in Figure 25. Both errors in- 
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crease as the actual range increases. However, the maximum error in temperature and 
range is insignificantly small. The maximum error in ranging is 0.09 % and the maximum 
error in temperature is 0.0008 %. This results shows that the method performs remark- 
ably well for the conditions tested. 

In the calculation of the transmission factor for the i-th band, the precise t, should 


be calculated according to the weighted average 


A upper 


; OA, T)tiy(A)a2 
= (5.27) 


Taper 
0,4 ’ T)dd 


Riese 


rather than the standard average (5.12). It should be noted that, the flatter the photon 
emittance distribution within the spectral bands, the smaller the discrepancy between 
calculations (5.12) and (5.27). 

Despite the simplifying assumptions used in the model, the results demonstrate the 
feasibility of extracting both target range and temperature from thermal radiation 


measured over three bands. 
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Figure 15. Definition of vertical line from the earth through the receiver 
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Figure 16. 
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Figure 17. Definition of layers in atmosphere 


a2 





Figure 18. Layer of the atmosphere 





Figure 19. Geometry of ranging 





Figure 20. Geometry of ranging (cosine law) 
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Figure 21. 
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Figure 22. 
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Figure 23. 
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Figure 24. 


Error of range determination 
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Figure 25. 
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VI. CONCLUSIONS 


Previous work on dual spectral band measurements of radiant emittance demon- 
strated that, under somewhat ideal conditions, the ratio of the measurements could be 
used to determine the temperature. The main stipulated condition was that the bands 
were not self-contained. In the Chapter III of this thesis, the same conclusions were 
obtained for measurements of photon emittance. 

A general model encompassing arbitrary wavelength-dependent emissivity, 
responsivity, and transmission factor was proposed in Chapter 1V. Although the general 
problem is not analytically tractable, the methods applicable to constant emissivity, 
responsivity, and transmission factor (Chapter III) are approximately valid if the spec- 
tral bands are narrow enough. However, from the practical view point of making a 
signal measurement in the presence of noise, the spectral bands should not be too nar- 
row. Therefore, the more general model is often needed under practical conditions. For 
reason explained in Chapter IV, the photon detector responsivity introduced 1s fairly 
insensitive to wavelength up to the cut-off wavelength. In addition, for a broad class 
of materials emissivity curves are quite flat. Therefore, the focus of the reminder of the 
thesis concentrated on effects due to atmospheric attenuation. Specifically, knowledge 
of the range to the target would be necessary to apply the method discussed in Chapter 
IV. 

The concept presented and tested in Chapter V is created from the requirement that 
the dual spectral band method of temperature determination would require range infor- 
mation to determine the temperature. This assumes a disparity in attenuation effects 
between bands. By employing three bands instead of two, this problem 1s circumvented. 


Additional complications created by changing atmospheric conditions with trajectory 
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path are handled within the model. The numerical atmospheric attenuation data used in 
the model testing was taken from LOWTRAN. The results for this preliminary model 
are extremely encouraging. Future work needs to be performed to develop and test a 


more rigorous model. 
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APPENDIX A. DERIVATIVE OF & WITH RESPECT TO T 


The purpose of this appendix is to prove equation (3.7) which appears in the main 
text. Equation (3.7) had been introduced in [Ref. 8: pp. 1257] without proof: Since the 
ratio of photon emittance is a function of a, it requires a few steps to get its derivative 
with respect T. For the reader’s convenience, relations (3.2a), (3.5) and (3.6) are re- 


produced, which show that 


[wav 
a, (Al) 
[“wau 
where 
y? 
g(v) Seto (A2) 
and 
aa, (43) 


Following the product rule, the derivative of the ratio of photon emittance (AI) with 


respect to temperature is 
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won a2 ~ Us ar _O_ oa 
wa [sca] eran — [ownav [acne 


éT aX 2 
| oa 


where the temperature dependence 1s expressed through a. One approach to simplifying 


(A4) 


the expression is to transfer the dependence of the integral limit on a evident in (A4) to 


the integrand. A substitution, 
W=ax (A5a) 


will transfer this dependence. It follows that 


x; 


[Peway =a] “elaxae. (45b) 


Therefore according to the product rule, 


Z| eyyav = | “elaxidx +0 | F glax)dx. (46) 


aX 


Now from (A2) and the transformation rule (ASa), 


(ax)? 


exp(ax) — 1’ ba) 


g(ax) = 


and it follows from (A7) that 


ax? ax? exp(ax) 


exp(ax)—1 [ exp(ax)-—1]? ae) 


g = 
Ou g{ax) a 


Direct substitution of (A8) into (A6) leads to: 


a [% x; 7 ax)? ” (ax) exp(ax) 
an gay = | a ah «| Poor 


x; { 


Noting the equivalence between g(ax) defined in (A7), and the integrand of the second 


integral (A9), it follows that : 


ax x a 3 
-< [sway - 3] “elas — | eameaay dx. (A10) 


t 


After reversing the transformation (A5a), the integrands in (A10) will depend on w and 


(A10) leads to: 


a ae ae 3 
& | Pewnav = 3 [ewan - | toma 


For compactness of representation, define 


2 
fW)= yw exp(y) 


=O, Al2 
[ exp() — 1)’ ai 


65 


Then combining the result (A3) with (A11) leads to : 


oF | away = [Away -3 [away | (A13) 


After substituting (A13) into (A4) 


[nay —3 [Pew [ecnav - [ewe] [nna -3 Per 








0B, I 
ao oP aX, ; 
[“ewev 
aX; 
(A14) 
which after simplification leads directly to (A16) or (3.7), as given in the text. 
aX AX, aX, oN 
["nvrav[ewrav — [eav | “nna 
02, ] a aX; a aX3 
= (A15) 


| “sca 


aX; 
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APPENDIX B. ASYMPTOTIC APPROXIMATION FOR RATIO 
CALCULATION 


The purpose of this appendix is to prove equations (3.22) and (3.23) which appear 


in the main text where the asymptotic approximations are applied. 


A. FOR LARGE ALPHA 
For large a (i.e., small temperature), since x; values are fixed, » takes on large posi- 


tive values Therefore, g(y) and fw) functions are approximated from (3.5) and (3.8) as 


follows: 
2 
Wi soy (BI) 
and 
y? 
IDES Seay « (B2) 
After direct substitution into (3.8), it follows that: 
ax? aX, 
va Wy exp(—Waw|  w’ exp(— Wd 
. a ee | (B3) 
- | Wy exp(—W)dy| Ww exp(— pd 


One round of integration by parts yields 
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aX aX, 
veil — W exp(— ¥) J? 43 | y’ exp( — vv} | Wy exp( — Wd 


” o a (B4) 
- | y’ exp( — vavt —W exp(—w)]R+3) yw’ exp(— ya} 
After algebraic simplification, (B4) becomes 
BX 
ifgl=[ — ¥° exp(— ni] py” exp( — Wd 
% (B5) 


-| Wy’ exp(— Wdb[ — W° exp(— v) ox 


4 


According to [Ref. 12: pp. 85] 


| y’ exp( — Wd =— w’ exp( — p) — 2 exp( — #) — 2 exp(— y), (B6) 


which, for large positive w (B6), is approximately given by 


\v exp( — W)dp~ — W’ exp(— y). (B7) 


After substitution of the approximate form (B7) into the integrals (B5), [f,g] becomes 
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igl= Ly’ exp( — WI" Lv? exp( — W188 — Cy? exp( — WIR"LW? exp( — ¥) ES. (B8) 
which, after simplification, becomes 


[fgl=L(ax,)° exp( — ax.) — «° exp( — a) JL (ox4)* exp( — x4) — (ax3)’ exp( — ax;)] 


— [(ax,)? exp( — x.) — a” exp( — a) IL(ax,)* exp( — ax4) — (ax3)° exp( — ax;) 1 


Since, for large a, the second terms within each bracket of (B9) become dominant, [f,g] 


is approximated as 
[Agl=o" [x3 exp( — a.x3) exp( — a) —x} exp( — a.x;) exp( — «)]. (B10) 
Recombination of (B11) leads to: 
[fg}xo°x3(1 — x3) exp [ — a(x; +1)}. (B11) 


This last form is applied in Section C of Chapter III. 


B. FOR SMALL ALPHA 
For small « (1.e., large temperature ), ¥ is small. Thus, g(yW) and f/(W) are approxi- 
mated from (3.5) and (3.8) by binomial expansion. After application of binomial expan- 


sion, g(wW) is integrated as 


| g(W)dy~ | u(1 + ay. (B12) 


And a binomial approximation for small w leads to: 


[ewrav=[v(1-> av, (B13) 
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which 1s integrated as 


fewuvav'{ 4-5 ) 


(B14) 


Similarly, after applying binomial expansion to (3.8), the integration of f(w) becomes 


[avrav~[va + w(1--} Pav. 


Algebraic expansion of (B15) leads to: 


which is approximated as 


[anaaev'(4-4 v). 


Therefore, for small a, [f,g] becomes 


“ nae ica) 


x 


which 1s expanded as 
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(B15) 


(B16) 


(B17) 


(B18) 


5 7 
3 
isl py Wa Sy) = 5 (> —1) (4 — 3) — 5 (2 —1)a — 23) 
2 dy 3 SiG aed a! 3 eh) 
: 4 4 
— Ty — x3 )(x3 —1) +35 ° (x4 —%3 (xy -1)- 37% = %3 xy = 1) 
Since, for small alpha, o5 terms are dominant, [{f,g] approximated as 
a ag Re 2 3. 3 
ifg]= i L(x eC) Call) i (x4 — x3) J. (820) 


This final form is apphed in Section C of Chapter III. 


71 


APPENDIX C. ALGORITHM TO CALCULATE THE RATIO OF 


PHOTON EMITTANCE 


band limits 1,2,3,4 


temperature range T1,T2 
criterion for error CR 













T=7T1+naT 





photon emittance 
calculations 

band limits 1,2 
band limits 3,4 






yes 


ratio calculation 
<> 


es 


integration by photon or radiant 
Simpson's rule orients Sal 


No 
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APPENDIX D. ALGORITHM TO CALCULATE THE RATIO OF 
OUTPUT VOLTAGE AT THE DETECTOR 


band mits 1,2,3,4 


















photon or radiant 
temperature range T1,T2 emittance distrioution 
criterion for error CR 
atmospheric 
T=T14+naAT transmittance 
photon emittance integration by detector 
calculations Simpson's rule responsivity 
band limits 1,2 
band limits 3,4 
no emissivity 


yes 





APPENDIX E. MODEL FOR PASSIVE RANGE AND TEMPERATURE 
DETERMINATION 


A. MAIN ALGORITHM 


detector heig Si 
Measured ratio data haat cual ! LOWTRAN data 
iota , 
layer 


Approximate 
Attenuation 
Coefficients 


| Measrerato cet guess P T 


new guess 
R,T 


FaNQe calculation 


temperature Ideal ratio 
see EQ(5.7) 
see EQ(5.25) 


No error < criterion 





yes 


range 
temperature 


B. LOWTRAN DATA GENERATION 


Inputs 

Wavelength Range 
Distance Range 
Atmospheric Conditions 
Angle=90 degree(vertical) 






LOWTRAN 6 Program 
Transmission Data for 
Each Layer/Wavelength 


Calculate Intrinsic Attenuation 


for Each Layer/Wavelength 
see EQ(5.11) 





C. LAYER (SUBROUTINE) 


Inputs 
Angle(Line of Sight) 
Range Guess 


ntrinsic Attenuation Coefficient 
Layer Thickness 


Detector Location 
Earth Radius 







Calculate r; 
see EQs(5.17),(5.21),(5.22) 


Calculate Average Attenuation 
Coefficients 
see EQ(5.15) 


o END 
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D. MEASURED RATIO PREDICTOR 


Inputs 

Range R1,R2 
Temperature T 
Angle 


R=R1+nAR 


R=Re 


Layer Subroutine 


Correct Average 
Attenuation Coefficients 
Measured Ratio 
Calculation 
Subroutine 


yes M 
<> Store R(T) 
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